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Nuclear pulse signal needs to be transformed to a suitable pulse shape to remove noise and improve energy 
resolution of a nuclear spectrometry system. In this paper, a new digital Gaussian shaping method is proposed. 
According to Sallen-Key analog Gaussian shaping filter circuits, the system function of Sallen-Key analog 
Gaussian shaping filter is deduced on the basis of Kirchhoff laws. The system function of the digital Gaussian 
shaping filter based on bilinear transformation is deduced too. The expression of unit impulse response of the 
digital Gaussian shaping filter is obtained by inverse z-transform. The response of digital Gaussian shaping 
filter is deduced from convolution sum of the unit impulse response and the digital nuclear pulse signal. The 
simulation and experimental results show that the digital nuclear pulse has been transformed to a pulse with a 
pseudo-Gaussian, which confirms the feasibility of the new digital Gaussian pulse shaping algorithm based on 


bilinear transformation. 


Keywords: Digital shaping, Bilinear transformation, Convolution sum 


DOI: 10.13538/j.1001-8042/nst.26.010402 


I. INTRODUCTION 


In a nuclear spectrometry system, the peak amplitude of 
nuclear pulse from the detector is proportional to energy of 
the incident particle, so the particle energy can be determined 
by measuring height of the pulse signal. A major cause of in- 
accuracy in the energy measurement is electronic noise, bal- 
listic deficit, etc [1]. Shaping techniques are intended to im- 
prove the signal-to-noise ratio and extract height of the pulse 
signal accurately [2-12]. Gaussian pulse shaping, as a gen- 
eral shaping method in the nuclear spectrometry system, per- 
forms good in this regard [12, 13]. With the rapid devel- 
opment of digital integrated circuits, and their advantages of 
flexibility and stability, studies on the digital Gaussian shap- 
ing algorithm for nuclear pulse signal has become popular 
during the past decades [12-21]. Valentin T. Jordanov stud- 
ied digital Gaussian shaping algorithm of exponential signals 
[15]. Chen et al. made efforts in digital Gaussian shaping al- 
gorithm based on wavelet analysis [13, 14]. Zhou et al. made 
advances in digital Gaussian shaping algorithm based on nu- 
merical differentiation [19]. The digital nuclear pulse signals 
are treated as quasi-Gaussian signals, and recursive algorithm 
is used to achieve the digital Gaussian shaping. In this paper, 
a digital Gaussian shaping algorithm based on bilinear trans- 
form is proposed. The unit impulse response of the digital 
Gaussian shaping system is deduced, and Gaussian shaping 
is realized by calculating the convolution sum of unit impulse 
response and digital nuclear pulse signal. 
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II. SYSTEM FUNCTION OF THE ANALOG GAUSSIAN 
SHAPING SYSTEM 


The analog Gaussian shaping in the nuclear spectrometry 
system, which is often implemented in Sallen-Key filter, as 
shown in Fig. 1, can shape a nuclear pulse signal as a pseudo- 
Gaussian signal. 


Fig. 1. Scheme of the Sallen-Key filter circuits. 


According to Kirchhoff's law, the Eq. (1) can be obtained, 
[f (t) — y1(0]/R = Cdy»(0/dt + Cd O - y(0)/dr, — (1) 


| py) 
yx(t) = ROR (2) 
yi(f) = yo(t) + RCdy2(t)/dt. (3) 


Thus, the relationship between input signal f(t) and output 
signal y(t) can be expressed as 


RCY” E) + RCy (0) + y(t) = 2f (0). (4) 
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Let the Laplace transform of f(f) and y(t) be F(s) and Y(s), 
respectively. Using differential properties of Laplace trans- 
form, we have the Laplace transform of y'(r) and y’’(f) be 
sY(s) and s?Y(s), respectively, and we have 


R?C? s Y(s) + RCsY(s) + Y(s) = 2F(s). (5) 


Thus, the system function of the analog Gaussian shaping 
system is 


H(s) = Y(s)/F(s) = 2/2 C? ? + RCs + 1). (6) 


Il. SYSTEM FUNCTION OF THE DIGITAL GAUSSIAN 
SHAPING SYSTEM 


$,-plane z-plane 


Fig. 2. Mapping relationship of bilinear transformation. 


In the process of transforming Gaussian shaping system 


ChinaX ive ERAT 


Nucl. Sci. Tech. 26, 010402 (2015) 


from analog domain to digital domain, bilinear transforma- 
tion is used for conversion in order to maintain the basic 
properties of the system function unchanged and overcome 
the spectrum aliasing phenomenon. The principle of bilin- 
ear transformation is secondary mapping method, as shown 
in Fig. 2. 

The mapping process includes three steps, 

(1) The s-plane is mapped onto the strip-like region from 
-1/T to n/T of the s;-plane. A conversion can be done as 


s = 2tan(s,T/2)/T = (1 e *9T)/(1 e? )[2/T]. (7) 


(2) Mapping the s;-plane onto the z-plane; the conversion re- 
lationship can be expressed as 


qeu, (8) 


(3) Mapping the s-plane onto the z-plane; the following 
transform equation between s and z domain based on the 
Eqs. (7) and (8) can be obtained, 


s-ü-zb/ü-szby/T]. (9) 


This allows the one-to-one mapping relationship from the 
s-plane to the z-plane. Thereby spectrum aliasing phe- 
nomenon is eliminated. The system function of the digital 
Gaussian shaping system can be obtained by Eqs. (6) and (9) 
when letting M = RC/T 


H(z) = (2 +47! «27 ?)J/[AM? + 2M + 1 (2 8MP)z ! + (AM? - 2M  1z?]. (10) 


Leta = AM? 32M +1, b = 2- 8M?, and c = AM? - 2M +1, 


the inverse z transform is carried out to Eq. (10). Then the unit 
impulse response of the digital Gaussian shaping system is 


h(n) =(2k; /a){[-b — (P? — 4ac) ?]/Qa)" u(n) + (2k2/a){[-b + (P? — 4ac)! ?]/Qa)" un) 


+ (2(k3/a){[-b — (b? — 4ac) ?]/Qa)"  u(n — 1) + ([2ka/a){[-b + (b? — 4ac) ?]/(2a)) ! u(n — 1), 


where, kı = (b — 4a)/[2(b* — 4ac)!?] + 1/2, kn = (=b + 
4a)/[2(D? — 4ac)!?] + 1/2, k = -a/(b* — Aac)!?, and k4 = 
a/(b? — 4ac)!?. 

Let z = e/^, the discrete-time Fourier transform of h(n) can 
be obtained by Eq. (10), 


H(e/") = (2+ 4e + 2e ?/")/(a + be + ce). (12) 


The amplitude spectrum of H(e/^ ) is shown in Fig. 3. The 
spectral width of the digital Gaussian shaping system increas- 
es with M decreasing. This makes the high-frequency noise 
filtration reduce. 


(1) 


IV. GAUSSIAN SHAPING OF THE NUCLEAR PULSE 
SIGNAL 


A. Mathematical description of the nuclear pulse signal 


The structure of a nuclear spectrometer system with digital 
shaping is shown in Fig. 4. If the charge generated by an im- 
pinging particle is instantaneously collected by the detector, 
the output signal from the preamplifier can be approximately 
expressed by a continuous-time exponential signal 


f(t) = Bee!" — Bye"! u(t), (13) 
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Fig. 3. The amplitude spectrum of H(e/”). 
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Fig. 4. Scheme of a spectroscopic system with digital shaping sys- 
tem. 


Detector +> 


where, B is magnitude of the signal; u(t) is a unit step signal; 
Tı and T2 are time constants of the exponential signal; and 
B, is proportionality coefficient. To simplify the calculation, 
let Bj = 1. If the ADC sampling period is T, the sampled 
discrete-time signal f(n) is 


fn) = B(e "7/1 — Bye ™ u(n). (14) 


B. Digital Gaussian shaping of the simulated double 
exponential signal 


The response of the input signal f(n) after the digital 
Gaussian system is described as 


y(n) = f(n) * h(n) 


= >) fhn- m), ml 


where, h(n) is the unit impulse response of the system; and 
y(n) is the zero state response of the system. 


Take double exponential Gaussian signal as input of the 
digital Gaussian shaping system, and the sampling period 
T = 0.01 us, B = 1.6, t; = 0.23 us. Fig. 5 shows the re- 
sponse y(n) of digital Gaussian shaping system, calculated 
and plotted by Matlab at t) = 0.4 and 1.8 us for different 
of M values. Width of the shaping waveform increases with 
M, and the symmetry improves with larger M. An advantages 
of the digital Gaussian shaping system is that one needs just 
changing the M value of the parameter, rather than adjusting 
the hardware. 


Amplitude (V) 


Amplitude (V) 
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Fig. 5. (Color online) Response of the digital Gaussian shaping sys- 
tem with simulation exponential signals. 


Amplitude (V) 


Time (us) 


Fig. 6. (Color online) Response of the digital Gaussian shaping sys- 
tem with sampling signals from detecting a Co source. 


C. Digital Gaussian shaping of the real sampled nuclear pulse 
signal 


Signals from a Co source were measured by a Nal(TI) 
detector and an Agilent DSO-X 2012A oscilloscope. The re- 
sponse y(n) of the digital Gaussian shaping system is calculat- 
ed by Eq. (15) in MATLAB, after the digital nuclear pulse sig- 
nals are stored, with sampling frequency of 200 MHz. Fig. 6 
shows the measured pulse signal and its corresponding wave- 
form of the digital shaping signal at M — 80 and 160. The 
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Fig. 7. (Color online) Response of the digital Gaussian shaping sys- 
tem with actual sampling signals at M = 0.5. 


shaping waveform at M — 160 is wider, and the symmetry 
is better, than those at M = 80. The shaping waveform is 
preferably Gaussian waveform. Decreasing the M value, the 
shaping waveform will be narrowed in time domain. Also, a 
low shaping time reduces the noise filtration and resulted in 
a poor energy resolution. At M = 0.5 (Fig. 7), significant 
noises in shaping waveform can be seen. 
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V. DISCUSSION AND CONCLUSION 


Pulse shaping digitization is the key to digitization of nu- 
clear instruments. After deducing unit impulse response of 
the digital Gaussian shaping system by bilinear transforma- 
tion, the convolution sum of the digital nuclear pulse signal 
and unit impulse response of the derived digital Gaussian 
system can be carried out. The digital nuclear pulse signal 
is shaped as pseudo-Gaussian signal and the noises are re- 
moved. A higher M leads to better symmetry of the shaping 
waveform, wider pulse, and better effect of noise filtering, 
with flat top pulse and small ballistic deficit. But a wider 
pulse will increase the pulse pile-up. Pile-up distortion is a 
major drawback in nuclear spectrometry system at high count 
rate. Narrowing the width of the pulse can lead to the reduc- 
tion of the pulse pile-up distortion, but a low shaping time 
reduces the noise filtration and results in a poor energy reso- 
lution. In practical applications, the shaping parameter values 
can be changed neatly to meet different shaping needs with- 
out adjustment of the hardware. 


The study provides a new method for achieving the digital 
Gaussian shaping of the nuclear pulse signal. 
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